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Abstract. Eractional calculus is a generalization of classical theories of 
integration and differentiation to arbitrary order (i.e., real or complex 
numbers). In the last two decades, this new mathematical modeling ap¬ 
proach has been widely used to analyze a wide class of physical systems 
in various fields of science and engineering. In this paper, we describe 
an ongoing project which aims at formalizing the basic theories of frac¬ 
tional calculus in the HOL Light theorem prover. Mainly, we present the 
motivation and application of such formalization efforts, a roadmap to 
achieve our goals, current status of the project and future milestones. 
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1 Motivation and Background 

Physical and engineering systems are classified as continuous, discrete or hybrid 
depending upon the nature of underlying system parameters. The rich theories of 
mathematics provide the necessary tools to study the behaviour of such systems 
ranging from very small biological organisms to the modern Quantum mechani¬ 
cal phenomenons. Generally, differential equations [39] and difference equations 
[T^ are used to characterize the dynamics of these systems. Consequently, the 
concept of higher-order differentiation and integration are widely studied in di¬ 
verse disciplines of science and engineering. For example, it is well understood 
that the first derivative (^/(Q) and second derivative (^/(t)) of a function 
describe the rate of change and measure of concavity, respectively. However, we 
rarely think what if the order (n) of higher-order derivative {-^) becomes a real, 
complex or an irrational number? One immediate question arises in our minds 
is the existence or possibility of such a concept in mathematics. Interestingly, 
this seemingly new concept dates back to 1695 when L’Hopital asked Leibniz 
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regarding his notation “what if n is In reply, Leibniz [20] prophesied 
in his letter, .. Thus it follows that d^x will be equal to xyfdxTx. This is an 
apparent paradox from which, one day, useful consequences can be drawn ... ”. 
Leibniz’s initial work on the problem of defining the derivative of arbitrary order 
gave birth to a new field of research in mathematics (called fractional calculus) 
and attracted the attention of many physicists, engineers and geometers. Some 
of the great mathematicians and physicists who touched the field of fractional 
calculus are Riemann, Liouville, Laurent, Heaviside and Riesz [24]. 

The concept of fractional calculus has great potential to change the way we 
model and analyze the systems. It provides good opportunity to scientists and 
engineers for revisiting the origins. We briefly outline some of the the main ap¬ 
plications of fractional calculus in Table The importance of fractional calculus 
can be realized by the following quote from Miller and Ross [24] . They stated: 

“... The fractional calculus finds use in many fields of science and en¬ 
gineering, including fluid flow, rheology, diffusive transport akin to 
diffusion, electrical networks, electromagnetic theory, and probabil¬ 
ity. ... It seems that hardly a field of science or engineering has re¬ 
mained untouched by this topic ...” 


Field 

Applications 

Control Engineering 

- System identification [17] 

- Biomimetic (bionics) control [7] 

- Trajectory control [11] 

- Temperature control [30] 

- Fractional controller [53] 

Signal Processing 

- Fractional order integrator [19] 

- Fractional order FIR differentiator |37| 

- IIR-type fractional order differentiator [38] 

- Modeling of speech signals [18] 

Image Processing 

- Image restoration and edge detection [29] 

- Satellite image classification [6] 

Electromagnetics 

- Fractional curl ooerators I13I25I 

- Fractional Rectangular waveguides [14] 

Communication 

- Secure chaotic communication [I] 

- Informational network traffic modeling [40] 

Biology 

- Neuron modeling [5] 

- Biophysical processes [8] 

- Modeling of complex dynamics of tissues [55] 


Table 1: Applications of Fractional Calculus 


Nowadays engineering systems exhibiting fractional order dynamics are in¬ 
creasingly used in some safety-critical applications such as control systems, sig¬ 
nal processing, electromagnetics and electrical networks (as listed in Table [^. 















For example, fractional meta-materials based devices are used to build sensitive 
military and defence equipments and electromagnetic stealth technology |21| . 
Considering these facts, it is quite interesting and important to build a logical 
reasoning framework which can be used to formally verify such sophisticated 
applications within the sound core of a proof assistant. In fact, proof assistants 
have been successfully used to formalize and verify some challenging and para¬ 
doxical mathematical results, e.g., the formal proofs of the Kepler Conjecture 
(Flyspeck project) [TH] and the Odd Order Theorem m_ 

In this paper, we present details of an ongoing projeciQto develop a formal 
reasoning support for fractional calculus in higher-order-logic theorem prover. 
This project was originally started at the System Analysis and Verihcation 
(SAVe) lalj^in 2010. Earlier formalization was done in the HOL4 theorem prover 
with the main focus on fractional operators for real-valued functions and the 
verihcation of fractional order electrical components. Later on, the scope of the 
project was expanded to formalize fractional calculus involving complex-valued 
functions due to its various engineering applications (as listed in Table [^. Cur¬ 
rently, we are using the HOL Light theorem prover due to the availability of 
rich multivariate analysis libraries including Harrison’s recent formalization of 
complex-valued Gamma functior0as well as the interesting related projects like 
Flyspeck m and the formalization of optics theories (i.e., ray, wave, electro¬ 
magnetic and quantum) |2]. 

The rest of the paper is organized as follows: In Section we briefly review 
some commonly used notations and dehnitions of fractional order operators. 
We provide an outline of the proposed formalization framework in Section 
Consequently, the current status of the formalization and future milestones are 
discussed in Section]^ Finally, we conclude the paper in Section 

2 Mathematical Framework of Fractional Calculus 

There are different notations available for fractional derivatives and integrals. 
We use J^f{x) and D"" f{x) for fractional integral and fractional derivative, re¬ 
spectively. In these notations, v is the order of integration or differentiation and 
a is the lower limit of integration. 

For every function (/ : C —>■ C); and for every number e K or C, J" and 
D” should be related to / by the following criteria [TO] . 

1. If f{x) is an analytic function, then J^f{x) and D""f[x) must also be an 
analytic function of the variable x and of the order v of integration or dif¬ 
ferentiation. 

2. The operations J^f{x) and D^f{x) must produce the same result as ordinary 
integration/differentiation when u is a positive integer. 
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( 1 ) 


3. The fractional operators must be linear. 

Ja[o:f{x) + Pg{x)] = aJ^f{x) + pj^gix) 

D'^ [a fix) + Pgix)] = aD'^ fix) + pO'^gix) (2) 

4. The operation of order zero must leave the function unchanged. 

J°f = f and D°f = f (3) 

5. The law of exponents must hold for integration and differentiation of arbi¬ 
trary order under sufficient conditions on function /. 

= and D^iD-f)=D^+^f (4) 

Fractional integrals and fractional derivatives are also referred to as Differinte- 
grals m and there are more than ten well-known definitions for Differintegrals 
[H]- We describe here two of them, which are most widely used in analyzing 
real-world problems. These are the Riemann-Liouville and Griinwald-Letnikov 
definitions, which are also equivalent for a wide class of functions m- 

Riemann-Liouville (RL) Definition: 

JUix) = “ t)'’~^fit)dt ( 5 ) 

where J^/(a:) represents fractional integration with order v and lower integration 
limit a. The parameter a = 0 gives the Riemann definition and a = —oo gives 
the Liouville definition of fractional integration. Indeed Equation ([^ is the gen¬ 
eralization of Cauchy’s repeated integration formula to non-integer v [32]. Where 
r (.) in the above definition denotes the Gamma function which is defined using 
the well-known improper integral as follows: 

pOO 

riz) = / e-^e-*dt ( 6 ) 

Jo 

for Reiz) > 0. 

The fractional differentiation is given as follows: 

D'^fix) = ( 7 ) 

where m represents the ceiling of v, i.e., |"u]. 

Griinwald-Letnikov (GL) Definition: 

[X^] 

cDlfix) = lim h"” ^ (-1)'= (^\ fix - kh) (8) 

k=0 ^ J 

Griinwald-Letnikov definition caters for both fractional differentiation and inte¬ 
gration, as positive values of v give fractional differentiation and negative values 
of V give fractional integration. Here, (^) represents the binomial coefficients, 
which are described in terms of the Gamma function. 



3 Formal Analysis Framework 


The proposed framework, given in Figurej^ outlines the main ideas and roadmap 
to formalize the basic theory behind fractional calculus. The whole framework 
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Fig. 1. Formalization Framework for Fractional Calculus 


can be decomposed mainly into three major parts which are the formalization 
of the core definitions of fractional order operators, formalization of supporting 
transformations (i.e., Laplace transform |26j . Z-transform |28] and R-Transform 
[S]) and engineering applications. The first part heavily relies upon the Gamma 
function as mentioned in Section So the core step is to formalize the Gamma 
function in higher-order logic (HOL) and verify its important properties. Gonse- 
quently, any definition of fractional order operators can be formalized in HOL. 
However, our focus is two main definitions, i.e., Riemann-Liouville (RL) and 
fractional difference which indeed represent continuous and discrete versions of 
fractional order operators, respectively. This step also involves the validation 
of all the properties mentioned in Section This requires some important re¬ 
sults of multivariate calculus such as the notion of Lebesgue measurability and 
Fubini’s theorem which provides the reasoning support for iterated and double 
integrals. Interestingly, both of these requirements are available in the multi¬ 
variate analysis libraries of HOL Light. The second step is the formalization 
of important integral transforms which are necessary to analytically solve lin¬ 
ear fractional differential and difference equations. We mainly focus on three 
transforms, namely the Laplace transform, the Z-transform, and the recently 
introduced R-transform [1]. All of these transformations are used to transform 
complicated fractional differential (or difference) equations to algebraic equations 
which are easier to manipulate and to deduce interesting properties. Building 












































































upon these fundamentals, our ultimate goal is to formally verify a variety of en¬ 
gineering systems including control systems, signal processing, electromagnetics 
and electrical networks. All formalization steps make use of different multivariate 
theories of HOL Light, e.g., derivatives, integrals, complex vectors and measure 
spaces. Finally, the developed libraries of this project will become part of the 
existing HOL Light libraries. 

4 Current Status and Future Milestones 

As mentioned earlier, the project initially considered only real-order fractional 
operators in the HOL4 theorem prover. The main difficulty was the little support 
to handle improper integrals which was required to formalize the real-valued 
Gamma function. Therefore, we extended the integration theory of HOL4 by 
formalizing a variant of improper integrals using sequential limits. This was 
then used to formalize the Gamma function and verify some of its main proper¬ 
ties, such as the pseudo-recurrence relation [r{z + 1) = zr{z)), the functional 
equation (T(l) = 1) and the factorial generalization {r{k -I- 1) = k\) [34j. We 
utilized these foundations to formalize Differintegrals, given in Equations (H 
and 0. which in turn can be used to represent the dynamics of fractional or¬ 
der systems in higher-order logic. We also verified theorems corresponding to 
some commonly used properties of Differintegrals namely Identity and Linearity. 
Gonsequently, we conducted the formal analysis of a fractional order electrical 
component namely resistoductor, a fractional integrator and a fractional differ¬ 
entiator circuit [33| . Later on, the scope of the project was revised to include 
complex-valued functions and complex order fractional operators in HOL Light. 
The main requirement was to formalize the complex-valued Gamma function, 
Laplace and Z-transforms. However, the Gamma functiorj^ was formalized by 
Harrison in early 2014. In the meantime, we formalized the basic theories of the 
Laplace transform |36j and the Z-transform |35j . Currently, we are working on 
three main topics which include: 1) formal proofs of the uniqueness of Laplace 
and Z-transforms which are required to formally verify the inverses of these 
transforms; 2) vectorial Z-transform, which extend the simple Z-transform over 
complex vectors; and 3) fractional difference equations, which are mainly based 
on Gamma function, infinite summations and products over complex functions. 
Finally, we outline the major tasks to achieve the future milestones as follows: 

— Formalization of R-Transform. 

— Formalization of Differintegrals for complex-valued functions. This is mainly 
the generalization of the formalization which was developed in HOL4. 

— Formalization of linear fractional differential and difference equations with 
support to analytical solutions using the transform methods. 

During the course of this project, two master and two PhD students have 
contributed to the formalization. Interestingly, all of them are mainly electrical 

^ https://code.google.com/p/hol-light/source/browse/trunk/Multivariate/ 
gamma.ml 




engineers without prior background of formal methods and higher-order-logic 
theorem proving. Given the complexity and interdisciplinarity of this research 
project, it is quite encouraging to see people with an engineering (or physics) 
background to use proof assistants as a complementary tool. The formalization 
of the fractional calculus is quite challenging as it requires advanced mathe¬ 
matical concepts of vector integration and Lebesgue measurable functions, etc. 
So expertise in formal reasoning about these complex mathematical phenomena 
is required for this formalization, which is quite unique compared to reason¬ 
ing about software and digital hardware systems. The learning curve of HOT 
Light varies from student to student. Generally, students start proving basic 
math equations after a couple of months and the pace of formalization increases 
over time. Learning HOL Light libraries is not difficult once the basic concepts 
have been grasped by the user. The formalization of the improper integrals, 
the Gamma function, the fractional calculus, the Z-transform and the Laplace 
transform is approximately 15,000 lines of HOL Light code. One of the major 
obstacles in the formalization was the identification of suitable mathematical 
definitions and models. Sometimes, textbook proofs do not follow due to various 
reasons (corner cases, or the proof steps are too abstract, etc.) and they needed 
to be re-proved on paper with subtle details. Gonsequently, we have to modify 
the definitions and thus change the proofs. But now the current formalization 
seems quite stable as most of the classical properties have been formally ver¬ 
ified for our definitions. We believe that future developments can be built on 
the foundations that have been formalized as most of the work is for general 
systems. Finally, another important aspect of this project is the potential to 
apply developed theories to various applications other than fractional calculus. 
For example, we demonstrated the use of the Gamma function in probability 
theory [34] , the Z-transform in signal processing [35] , and the Laplace transform 
in power electronics [36]. 

5 Conclusion 

In this paper, we mainly presented the motivation and ongoing activities of our 
long term project about the formalization of fractional calculus in the HOL 
Light theorem prover. The main contribution of this project is a comprehensive 
framework of formal definitions and theorems about fractional calculus which 
can be used to verify modern control, signal processing and electromagnetic 
systems. Some future directions and recommendations for HOL Light are the 
improvements in the visualization of proofs, better automation and more acces¬ 
sible tutorials with examples from different engineering/physics topics. 
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